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We present a conjecture describing new long range correlations among the Riemann zeros 
leading to 3 principal features: (i) The spectral auto-correlation is invariant w.r.t. the aver- 
aging window, (ii) Resurgence occurs wherein the lowest zeros appear in all auto-correlations, 
(iii) Suitably defined correlations lead to predictions [prophecy) of new zeros. This conjec- 
ture is supported by analytical arguments and confirmed by numerical calculations using 
10 22 zeros computed by Odlyzko. The results lead to a self-duality of the Riemann spectrum 
similar to the quantum-classical duality observed in billiards. 
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PACS numbers: 02.10.De,03.65.Sq, 05.45.Mt, 05.45.Ac 

I. INTRODUCTION 

The Riemann zeta function is 
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n=l p 

with product over all prime numbers p. The Riemann zeta function has a simple pole at s = 1. 
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Besides the trivial zeros at s = —2, —4, —6, • • • , the Riemann hypothesis (RH) states that all the 
nontrivial zeros are distributed on line as | ± itn with > hi. The Riemann zeta function and 
its zeros have been studied for over a century 0, 0, H 0, 0, 0,M II ED, E3, E1JI1 El, EI 13, 03 • 
The RH still remains unproven and presents a major mathematical challenge |l8j |. 

Motivated by the Hilbert-Polya conjecture of the existence of a Hamiltonian whose quantum 
spectrum is comprised of the set {t n } as the spectra, there has been a continuing quest to find 
such a Hamiltonian since it would prove the RH. Attempts have been made to construct certain 
■ Hamiltonian such as that with a fractal potential llflll for the Riemann zeros (RZ). For a brief 

^ , review about Hamiltonians and prime numbers, see |2fJ. 

The RZ {\±it n } can also be considered as a model spectrum for arithmetic or quantum chaos 
21 1 . It is widely believed that the ordinates {t n } of the RZ can be interpreted as the quantum 
spectrum of a chaotic system which does not have time reversal symmetry 0, though other 
interpretations are also possible. The most notable result in support of this notion is that of 
Montgomery who showed that the spectral correlation function R2{s) of the RZ is identical 
with that of a Gaussian Unitary Ensemble (GUE) which describes a chaotic system with no time- 
reversal symmetry in random matrix theory. This result has led to a line of inquiry which studies 
the statistics and correlations among the RZ. The merit of this approach is that the properties of 
the RZ will lead to insights into other systems where the classical dynamics is well established. 
These are motivated by the analog of the Riemann zeta function with Gutzwiller's trace formula 
and the dynamic zeta function of chaotic systems. 

In this paper, we present a conjecture which reveals new long range correlations among the 
RZ, obtained from a study of the auto-correlation of a smoothed Riemann spectrum. Several new 
properties emerge: (i) The spectral auto-correlation is invariant w.r.t. the averaging window. We 
demonstrate numerically that this remarkable invariance applies over any part of the spectrum of 
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10 zeros, and clearly suggests that the same structure is encoded in all parts of the spectrum, 
(ii) Resurgence occurs wherein the lowest zeros appear in all auto-correlations. This property of 
resurgence was first pointed out by Berry and Keating (hi) Suitably defined correlations lead 
to predictions (prophecy) of new zeros. 

The property of resurgence and prophecy suggests a self-duality for the RZ which we represent 
as C{\ ± it n } =>■ {\ =t tin}- This is very similar to a duality C{k n + ik' n } {7^ + irf} which is 
observed from the auto-correlation of the quantum spectrum of open chaotic billiards, and which 
relates the quantum resonance s pec trum {k n + ik' n } and classical resonance spectrum {7^ + ij"} 
of the hyperbolic n-disk repeller [2J|. 

The paper is organized in the following. In Sec. II, we introduce a smoothed spectral function 
and the spectral correlations and describe a conjecture based upon the spectral correlations. A 
"semiclassical" derivation of the conjecture is presented in Sec. III. Numerical calculations support- 
ing the properties of invariance, resurgence, prophecy and self-duality embodied in the conjecture 
are presented in Sec. IV. In Sec. V, we discuss the self-duality of the RZ with the quantum- 
classical duality observed in several billiard systems, and discuss the implications for a dynamic 
interpretation of the RZ. 



II. CORRELATIONS AND A CONJECTURE 

The spectral density is the sum of ^-function p(k) = ^2 n 8(k ± t n ). In order to study the 
correlation among the RZ, we use a Lorentzian-smoothed spectral density 



tt ^ (k ± t n ) 2 + e 2 ' 

Here t n are the ordinates of the RZ and e > is a small width. In the limit e —* 0, one gets the 
stick spectrum. It can be written as a continuous part plus fluctuations p e (k) = (p e (k)) + Sp e (k). 
The continuous part of the spectral density for the nontrivial RZ [J] is (p(k)) = (1/2tv) \n(k/2ir) in 
the limit e — > 0. We define the following fluctuation part of the spectral density |2j| 

Sp e (k) = -ljL[mn((l+e + ik)} 
tt ak 2 

2tt 2tt 7r ^ 

n=l 



+ 



( fc + tn )2 + e 2 {k _ tn) 2 +e 2 



(1) 



The function 5p e (k) defined above is an even function. Note that the above quantity Sp e (k) is related 
to the fluctuation part of the quantum time delay tq(w) considered in [2^, Spi/ 2 (k) = iTTf[(k/2)/A. 
In practice, 5p t (k) is a finite sum. For the spectral function 5p e (k) in a window k £ [Kq — A, Ko + A] 
with Kq S> A 3> 1, one notices that the contribution of RZ far outside [Kq — A, Kq + A] is negligible 
for small e. Thus one may use a finite sum 

1 k 1 N e 

5pJk)~ In h- ^ 71 tt5 o- (2) 

H y ' 2tt 2tt tt ^ (k-t n ) 2 + e 2 y 1 

n=N y n> 

The RZ used in the above expression should be slightly outside the window such as tN % Kq — A 
and t N > K + A. 

Define the following cross correlation of 8p e {k) among two windows [K\ — A,Xi + A] and 
[K2 — A, K2 + A] with K\, if 2 any real numbers and K\ < K2 as 

C^( S ;^i,K 2 ,A) = {{Sp^ + tySp^ + k + s)}^, \s\<A, 

k G [-A — min(0, s), A — max(0, s)]}. (3) 
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The superscript "c" denotes cross. Here the average over k is defined as 



b — a 



h{k)f 2 {k)dk. 



There are two special cases 



C r pe (s; K , A) = C c pe (s; K , K , A) = (Sp e (K + k)8p e (K + k + s)) k , (4) 
C* (s; K , A) = C c pe (s; K , -K , A) = (5p e (K - k)Sp e (K + k + s)) k . (5) 

In Eq. (jSJ), we used Sp e (Ko — k) = 5p t (—Ko + k). The superscript "r" and "p" denotes resurgence 
and prophecy, respectively. The meaning will be apparent later. 

We now summarize the principal results of this paper. We have the following 

CONJECTURE: 

1. INVARIANCE the spectral correlation is invariant and independent of averaging window. 
That is 



C£(a;tfi,Jr 2 ,A) = C^-i^ + s), 
C; e ( S ;K ,A) = C e (s), 
C£>;if ,A) = C'(2K + s) 



(6) 
(7) 
(8) 



Note that one has —A < s < A with A (5 = t n — t n -\), the mean spacing between the RZ in a 
window. 

2. The spectral correlation is a sum of functions of the pole and zeros of the Riemann zeta 
function 



C £ (s) = R 



g(l + 2e + is, 1) - ^ 5(1 + 2e + is, -2n) - ^ g{l + 2e + is, - ± it n ) 



n=l 



n=l 
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Here <?(z, a) is a certain two-variable function. The second (third) term is the summation over 
trivial (nontrivial) RZ. 

3. The function g(z, a) can be broken into diagonal and off-diagonal terms 



g(z, a) = f/dia g (>, a) + tfoffO, «)• 
Here g Q s <C <?diag- Then we show that 

1 1 



)_2 



2vr 2 (z - a) 2 2vr 2 ^ (z - a/nf 

\ 1 n=2 v ' ' 



(10) 



(11) 



with c n 's real coefficients, \c n \ < 1/n, and c n — > for n — * 00. For the correlation function of 
interest C Pe (s) with s 6 [a, 6], it can be replaced with a finite sum over the pole and zeros of the 
Riemann zeta function 
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(is + 2e) 2 ^ (is + l+2e + 2n) 2 ^ (is - it n + ±+2e) 

n=l n=No * ' 

s G [a, b], t No <a, t N > b 

The contribution of RZ far outside [a, b] is small. 



(12) 
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FIG. 1: A sketch of the conjecture: (a) cross-correlation C r p ^(s), (b) invariance and resurgence in auto- 
correlation Cp e (s), (c) prophecy C?(s). 

4. RESURGENCE 



1 N 



^ (is ± it n + i+2e) : 



;, \s\ <A,t N >A. 



(13) 



5. PROPHECY 



1 - 

CP e (s;K , A) = C e (2if + «) « £ 



1 



n=N 



(is + 2iK - it n + 4+2e) 



2 > 



|s| < A, t No <2K -A, t N > 2K + A. 



(14) 



A graphic representation of the conjecture is shown in Fig^ 
A few remarks are in order. 

INVARIANCE The invariance of the correlations means that the same structure is encoded in 
the Riemann spectrum independent of the center and width of the averaging window. The cross 
correlation C c Pf (s ; K\ , if 2 , A) between the spectral function 5p e (k) in any two windows [if 1 — A , if 1 + 
A] and [if 2 — A, if 2 + A] is the same if the center distance if 2 — ifi is the same. 

RESURGENCE We find that the peaks of — C^(s) are located at the positions s = t n . The 
correlation between RZ in [ifo — A, ifo + A] leads to RZ in [—A, A]. The auto-correlation C pi (s) 
in Eq. (j3J) is essentially the two- level correlation i?2(s) 0. Resurgence was first noted by B erry 
and Keating 0, 0, 0] and has been revisited recently by Leboeuf and collaborators 0, 0, 0]- 
While the correlations studied in the present paper are related to the above work, the difference 
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is that here, there is no unfolding and no rescaling by the mean level spacing, and reveals some 
unexpected new results among the RZ summarized in the equations above. 

PROPHECY can lead to the prediction of new zeros in the interval [2Kq — A, 2Kq + A] from 
interval [K - A, K + A]. 

Resurgence and prophecy can also be observed in the cross correlation Cp e (s; K\, K2, A). 

III. TRACE FORMULA DERIVATION 

The above properties are "derived" using a semiclassical method 28]. The fluctuation part of 
the spectral density of the RZ is 0, Ei| 

1 00 1 

s Pe (k) = — Yl l p E a p r{l+t) c ° s ( rkL p) ( l5 ) 

p r=l 

Note that this is exactly of the same form as the Gutzwiller trace formula for chaotic billiards 



in wave number k 27j, with length of "periodic orbit" equal to the logarithm of prime numbers 



L n = Inp and A p = e Lp = p. Like that of the chaotic systems such as the chaotic n-disk systems 



the number of "periodic orbits" for the RZ grows asymptotically as N(L P < L) ~ e L /L. 
Although the above trace formula converges only for e > ^, we may still use it formally for e — > 0. 
The correlation can be written as the sum of diagonal and off-diagonal terms 



C;(s) = (5p e (k)5 Pe (k + s)) = Q iag (s) + C* s (s). 

The diagonal part is 

1 00 1 B 2 00 1 

<W*) = pE L ?E v (1+2e) ™< rsL P ) = - 2^*5? E E -2 f P ( 16 ) 

p r=l p r=l 

with t p = Ap (1+2e) e isL p = p- 1 -" 2 ^™. C| iag (s) can also be broken into two parts [HQ 

CkgOO = C e c (s) - CUs) (17) 



with 



CHs) = - J-R^ln^^). (18) 



Here 



h ^( s )=EE^; +1 ' (19) 

p r =l V ' 

Alternatively, one has the following expressions by direct summation over r 



Cdiag( s ) - 7^2^E L Pl_/ ' 

P P 

C W = i^E% tp 



(i-i P ) 2 ' 

<*M = ^»E%^W <») 
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These expressions have absolute convergence. Obviously, one has CJr(s) < jC^s). 

In order to obtain a compact form for C^(s), we notice the following expression for the Riemann 
zeta function £(s) 



e s/ln 



e (ln 2.-l- 7 /2) S / s X 

CW = 2( i -i)r(i + i a) n( 1 -^J 

with 7 n = 2 i ii n . One has 

£ ln c(1 + ls) = h + + E (i+u-^r (21) 

Here ip n ( z ) = ^ n+1 r(z)/dz ra+1 is the polygamma function with 

oo 1 

= (-1)™+^! V — -. 



Define the following function 



- OO ^ OO 



(22) 



one thus has 



q(s) = ^»Z(l + 2e + ta). (23) 



In order to obtain a compact form for the term C%(s), we use the following expansion for small 

t 

oo oo oo nr 

e^'=e^eV- (24) 

g=2 y n=2 r=l 

The coefficients c n can be obtained by equalizing the coefficients of t q on both sides of the equation. 
For prime numbers and the power of 2, one has 

c p = (p- l)/p 2 , 
c 2 n = 2" 2n , n > 1. (25) 

The rest of the c n 's are given by the following recursive equation 

E = ^ (26) 

n>2,r>l,nr=g 

For example, one has C6 = —1/18, c\q = —1/25, C12 = —1/72, c\4 = —3/98, C15 = — 8/15 2 . 
Using the expansion (|24|) and function Z(z), one has 

1 a2 00 

Cr(s) = -^»^£^ln<[7i(l + & + «)] 

n=2 

- 00 

= ^^^2n 2 c n Z[n(l + 2e + is)}. (27) 

n=2 
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One thus gets 

1 °° 

Cdiag(s) = -^^2n 2 c n Z[n(l + 2e + is)} 
n=l 

= -^2n 2 c n C^ +n \ns) (28) 



71=1 



with c\ = — 1. If one ignores the off-diagonal term Co ff (s), one obtains © with <?(z, a) given by 
Eq. (JTTJ). 

The value of the correlation at the origin can also be calculated in the following 



^ ag (o) = iE^Ev (1+2e) = iE- 



T 2 



2f 2 ^ P ^ P 27T 2 ^ e (l+2e)L p _ 1 ■ 

p r=l p 

The average density of periodic orbits for large L is p(L) = d(e L / L)dL = (L — l)e L /L 2 . Since one 
has 



/ L e (l+26)L_i ; L 4 e 2 

Simply set L = In 2, one has C| iag (0) « (1 + 2eln2)/87r 2 e 2 4 £ . For example for e = 0.5, one has 
C e Ai (0) ps 0.0429. In the limit e -> 0, one has C| (0) ps l/87r 2 e 2 . 



IV. NUMERICAL RESULTS 



In order to check the validity of the conjecture, we construct 5p e (k) according to Eq. (j2j) from 
a finite sequence of t n such that t^o ^ ifo — A and > i^o + A. Define 5at = 27r/ ln(tAr/27r) 
which is the local level spacing in the neighbourhood of tjy. We use e ~ 5n- One has 5at ~ 0.669, 
0.257, 0.141, 0.134 for N = 10 5 , 10 12 , 10 21 , 10 22 , respectively. The RZ we studied are obtained 
from Odlyzko [3(|. We find the following properties of spectral correlations of the RZ. 



A. Correlation invariance 



For the first 10 5 Riemann zeros, the fluctuation part of the spectral density 6p e (k) is constructed. 
We equally divide them into 5 segments : [K n - A, K n + A] with K n = (2n - 1)A + 10, A = 7490, 
and n = 1, 2, • • • 5. Note that t 10 s = 74920.8. We then evaluate the auto-correlation C r pe (s) Q) for 
each segment. We find that without any normalization, all 5 of them collapse to the same curve . 
As a comparison, we also calculated the auto-correlation for every 2 x 10 4 of the first 10 5 RZ. The 
correlations are also the same. The relative difference between them is less than 0.1%. 

In the same way, we consider four groups of RZ: {1, 10 5 }, {10 12 + 1, 10 12 + 10 4 }, {10 21 + 1, 10 21 + 
10 4 }, and {10 22 + 1, 10 22 + 10 4 }. For the spectral correlations of the above four segments of RZ 
with the same width e, they also all collapse to the same curve without any normalization as shown 
in FigHJ 

In general, one would expect that the spectral correlation will depend on the choice of Kq and A 
in Eq.Q). Even though the spectrum is more congested around t n ~ Kq for larger n, they all have 
almost the same spectral auto-correlation. We find that the main feature of the auto-correlation 
is only determined by the choice of e and independent of Kq and A. So no matter how high the 
RZ located on the critical line, they are all correlated in the same way as the low lying RZ. The 
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FIG. 2: Invariance of spectral correlation C r p (s; K ; A) for RZ with e = 0.25: {1,10 5 } (dashed), {10 12 + 
1, 10 12 + 10 4 } (dash-dotted), {10 21 + 1, 10 21 + 10 4 } (dotted), and {10 22 + 1, 10 22 + 10 4 } (solid). All curves 
for different ranges arc almost identical. 

dependence on Kq and A will show up for very small width e such that e < S and very narrow A. 
This is shown in Fig|3]where we used e = 0.1 < £ 10 22 ~ 0.134. 

We point out that the diagonal approximation is very accurate for the Riemann zeros. Numer- 
ically calculated C pe (s) are almost indistinguishable from the direct evaluation of Cdi ag (s) given in 
Eq. Q2U|) using the first 350,000 primes for s < 500. For example for e = 0.25, 0.5, one has the 
standard deviation a = (C 2 ^) 1 / 2 / '{C r2 ) 1 ' 2 ~ 0.022, 0.031, respectively for the spectral correlation 
of the RZ {t n } with n ~ 10 12 . 

B. Resurgence 

The invariance of the spectral correlations implies the existence of certain intrinsic structure in 
the spectral auto-correlation. By resurgence of RZ we mean the appearance of RZ {| ±it n } with 
small n from the correlation of RZ with large n. We find that the peaks of —C r p ^{s) are located at 
the positions s = t n . 

We then use our conjecture (|12jl to reveal the structure within the auto-correlation. The nu- 
merically calculated C r pe (s; Kq, A) can be approximated very well by C^(s) without any fitting 
parameters. This is shown in Fig. 0] for width e = 0.25. The difference C r Pt {s) — C^(s) is very 
small as shown in Fig. EK with a = ((C r pe - C^) 2 ) 1 / 2 / (Cf) 1 / 2 = 0.174,0.155, for e = 0.25,0.5, 
respectively of the RZ {t n } with n ~ 10 12 for < s < 350. Better approximations can be achieved 
by including more terms in C| iag (s). Define the difference of a better approximation according to 
Eq. (21, 
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FIG. 3: Same as in FigHwith e = 0.1. 

one has the standard deviation reduced to a = 0.093,0.072, for e = 0.25,0.5, respectively. We 
further define 

AC 3 (s) = AC 2 + 2C; i+ .J3s;K , A), (30) 

then the standard deviation a = 0.041,0.036 for e = 0.25,0.5, respectively. This is comparable 
with the off-diagonal term. The behaviors of C r Pt — C|, AC2, and AC3 for e = 0.25 are plotted in 
FigEI 

Resurgence can also be observed in the cross correlation C pe (s; K\, K 2 , A) in Eq. © for two 
windows [K\ — A,K\ + A] and [K 2 — A,K2 + A] with K\ ~ 0. Suppose we consider the set 
{t n } with 10 22 < n < 10 22 + 10 4 . We use K x = 0, K 2 = {t 10 2 2+1 + t W 2 2+W 4)/2, and choose the 
window size A = 620 ~ (t 10 22 +10 4 — t 10 22 +1 )/2 in @ so that K2 + A ~ t 10 22 +10 4. The spectral 
function 8p e (k) is constructed according to Eq. (^Q) in windows [—A, A] and [K 2 — A,K2 + A]. 
The center distance between these two windows is K2. The cross correlation Cp s (s; 0, K2, A) with 
s G [— A/2, A/2] and e = 0.25 is shown in FigEI One can see that C c Pt (s; 0, K2, A) can be well 
approximated by C^K 2 + s) with a = {{C c pe - C^) 2 ) 1 / 2 /(C 2 ) 1 / 2 ~ 0.199,0.171 for e = 0.25,0.5, 
respectively. To improve the fitting, one can use 

C c pt (s;0,K2,A)^Cl(K 2 + s)-C p ' i (2s;K 2 ,A). (31) 

Here C v Pt (s; K 2 , A) with -A < s < A is the cross correlation © between [—K2 — A, —K2 + A] and 
[K2 — A, K2 + A]. This time, the standard deviation is reduced to a = 0.118,0.079 for e = 0.25,0.5, 
respectively. Similar results are also obtained for K2 ~ £1012 and K2 ~ tio2i • 

We point out that the small s behavior of C pc (s) is not a Lorentzian, contrary to the claim in 
[2^ |. According to (|12|) . one has 

C;» ^ (4e 2 - s 2 )/2vr 2 (s 2 + 4e 2 ) 2 (32) 

for s < 5. So that C p ^(0) ~ l/(87r 2 e 2 ). This value is very close to that at the origin for different 
auto-correlations. The behavior of C r p<L (s) for small s is shown in Figd for different segments of 
RZ. 
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FIG. 4: Resurgence of RZ in C p (s; Kq, A) for e = 0.25. Here C r p {s; Ko; A) is the numerically calculated auto 
correlation in [Kq — A, Kq + A] with Kq = (t w i2 +1 + t w i2 +10 i)/2 and A w (tioi 2 +io 4 — iio 12 +i)/2- The 
dashed curve is C^(s). Vertical lines are located at t n and coincide with the valleys. 
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FIG. 5: Resurgence of RZ in C r p {s; Kq, A) for e = 0.25. Here C p (s; Kq, A) is the numerically calculated auto 
correlation in [Kq — A, Kq + A] with K o — A w iio 12 +i and Kq + A w tio 12 +io 4 - ( a ) Cp( s ; ^o, A) — C|(s), 
(b) AC 2 (s), (c) AC 3 (s). 
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FIG. 6: Resurgence of the RZ in the cross correlation C p (s; 0, K 2 , A)) between [—A, A] and [if 2 — A, if 2 + A] 
with K 2 = (^io 22 +i + ^io 22 +io 4 )/2 and A w (f 10 22 +10 4 — iio 22 +i)/2- The dashed curve is CUK2 + s). 



The ratio C p (s)/C p (0) can be very well fitted by 



2\2 



(33) 



with fitting parameter ~ e for small s. This is shown in Fig|HJ As we mentioned before, 5p e (k) is 
related to the fluctuation part of the time delay tq(w), so C r p ^{s) is related to the auto-correlation 
of the time delay c(e; w c , A) discussed in j2f| with e = 0.5 and s = 2s. Thus for small s, c(e; iu c , A) 
can be better fitted by Eq. (|32j) with /? = 0.46 as shown in Fig|SJo. 



C. Prophecy 

By prophecy of RZ we mean the appearance of RZ ± zi n } with large n from the correlation 
of RZ with small n. We consider the correlation Cp e (s; ifo, A) in Eq. © for RZ {1,10 4 }. The 
prophecy of new RZ from Cp e (s; if 2, A) with e = 0.25 is shown in Fig^U The agreement with 
Ea.p4|) is very good. 

Prophecy can also be observed in the cross correlation C pe (s; ifi, K2, A) in Eq. Q for two 
windows [K\ — A, K\ + A] and [if 2 — A, if 2 + A] with K\ + A ~ 0. As in the above subsection, we 
consider the set {t n } with 10 22 < n < 10 22 + 10 4 . We use ifi = -A, K 2 = (t l0 22 +1 + 1 10 2 2+10 4)/2, 
and choose the window size A = (^1022+104 — iio 22 +i)/2 — 50 = 620.8. 6p e (k) is calculated in 
windows [— 2A,0] and [if 2 — A,if2 + A]. The center distance between these two windows is 
if 2 + A £ 10 22 +10 4. The cross correlation C Pe (s; A, if2, A) is shown in Fig^J If s < 50, one has 
resurgence of RZ since s + A + if2 < tio22 + iQ4. If s > 50, one has prophecy of new RZ. Similar 
results are also obtained for if2 ~ iio 12 an d K% ~ iio2i • 

The above long range correlation can be used to calculate RZ from known ones. When applied 
to prophecy correlation Cp t (s), new RZ can be obtained. The Fourier transformation of C Pt {s) is 



C(x) 



C p Js)e lsx ds 



—x^ e itn 



X — "fX 



(34) 



12 




0.25 0.5 0.75 1 1.25 1.5 

S 




-0.02 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

S 



FIG. 7: Small s behavior of C r p (s) with (a) e = 0.25 and (b) e = 0.5 for RZ: {1,10 5 } (dashed), {10 12 + 
1, 10 12 + 10 4 } (dash-dotted), {10 21 + 1, 10 21 + 10 4 } (dotted), and {10 22 + 1, 10 22 + 10 4 } (solid). The bold 
solid line is C* c (s) (Eq.©). 
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FIG. 8: Small s behavior of C r {s)/C r p {Q) with (a) e = 0.25, = 0.24 and (b) e = 0.5 , (5 = 0.46 for RZ: 
{1, 10 5 } (dashed), {10 12 +1, 10 12 +10 4 } (dash-dotted), {10 21 +1, 10 21 + 10 4 } (dotted), and {10 22 +1, 10 22 +10 4 } 
(solid). Here (3 is the fitting parameter in Eq. 



13 




100 



150 



FIG. 9: Prophecy of the RZ in C» (s;K , A) for RZ {1, 10 4 } with e = 0.25, K = 4941.4 , and A = 4920. 
The dashed curve is CU2Kq + s). The vertical lines are at s — t n — 2Kq. 
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FIG. 10: Prophecy of the RZ in the cross correlation C c Pe (s; A, K 2 , A) between [-2A,0] and [K 2 - A, K 2 + A] 
for e = 0.25. (a) C c Pe (s; A, AT 2 , A), (b) (s; A, K 2 , A) - C|(s + A + K 2 ). Here C^(s + A + K 2 ) is calculated 
from Eq. with £ n only up to iV = 10 22 + 10 4 . 



with 7 = 1 + 2e. C(x) will have peaks at x = rL p . Since the above expression is a sum over 
different exponentially decaying modes, one may use the Prony algorithm to get t n . The above 
method provides a way to calculate or at least to estimate new RZ. 
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V. DISCUSSION: SELF-DUALITY OF THE RIEMANN SPECTRUM 

The motivation for this work has come from experimental observations in the microwave trans- 
mission of open n-disk billiards which led to the observation of classical Ruelle-Pollicott (RP) 
resonances in the auto-correlation of quantum spectra of hyperbolic n-disk open billiards 
The result established a new approach to quantum-classical correspondence by demonstrating a 
correspondence between the quantum and classical resonance spectra of an open chaotic system. 
Applying the same procedures developed there to the Riemann spectrum, we have arrived at the 
results described in this paper. There are other parallels which are summarized below. 



1. In n-disk open billiards mentioned above [241 ] . the experimental trace was fitted well by the 
sum of Lorentzians T(k) = ^2 n b n /[(k — k n ) 2 + k' 2 ]. Here {k n + ik' n } is the spectrum of 
eigen-wave vector resonances, which are eigenvalues of the Helmholtz (or Laplace-Beltrami) 
operator, with Dirichlet boundary conditions (V 2 + k 2 )tp = 0, ifi = on dD. We showed that 
the auto-correlation of the transmission C(n) = (T(k)T(k + k)) are well-described as a sum 
of Lorentzians C(k) = ^ ch/[(k — 7J 2 + j" 2 ]- Here r y i + i"f'(aie the classical RP resonances. 
The coefficients c, are related to the eigen-functions of the Perron-Frobenius operator. The 
Cj are all positive, whereas the coefficients are negative for the RZ. For hyperbolic sys- 
tems, Ruelle has related these eigenvalues to the time-evolution of classical correlations |3l| . 
The above result can be viewed in the context of dynamic Ruelle zeta functions which 
are relevant to the n-disk billiard. The semiclassical Ruelle zeta-function is Q(—ik) = 

Y\ p 1 — e t{kL p +w^p/2) J y/| A p | A^" 1 whose poles are the quantum resonances {k n + ik' n }, 

while the poles of the classical Ruelle zeta-function C/?( s ) = Tip ~~ e~ sLp /Ap^j are the 
classical resonances {7' + ij"} (H, |U The product is over all classical primitive periodic 
orbits without repetition. Here L p is the length of periodic orbits with the Maslov index \x v . 
A p is the bigger eigenvalue of the Monodromy stability matrix for billiard systems. Math- 
ematically the quantum-classical correspondence can be viewed as a duality between the 
poles of the semi-classical and classical dynamic zeta functions through the auto-correlation 
procedure, i.e. C{k n + ik' n } =>• {7^ + ij"}. 

2. For the two-disk billiard (i.e. n = 2), the semiclassical resonances in wave vector space are 
k n = [mr-H(l/2) In A]/(.R— 2a) with A = a — 1 + y a(a — 2) the eigenvalue of the instability 
matrix in the fundamental domain and a = R/a. Here n is odd for A\ representation, n even 
and n 7^ for A2 representation The classical RP resonances of the two-disk system 
are j n = [In A + inir]/(R — 2a). In this case there is a direct 1-to-l correspondence between 
quantum and classical resonances. 

3. In the cases of the rectangle billiard and the equal-lateral triangle billiard, the auto- 
correlation Cp e (s) are shown [35] to correspond to the Fourier transform of the trace of 
the classical evolution operator, i.e. 

00 00 

cm - 4r£* = E E ^ - rL A = E E cos{rsLp) - 

p r=l p p r=l p 

There again the property of correlation invariance and resurgence are clearly demonstrated. 

4. For a particle on a surface of constant negative curvature, a self-duality exists for the quantum 
momentum spectrum, {p n } and the classical spectrum = | db ip n } [if]]. The quantum 
eigen-energy are given by E n = \ + p 2 . 
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Motivated by these observations, we examine a dynamic interpretation of the Riemann zeta 
function. One has the following expression due to Riemann 

= r nr-^ 2 ^os(yn,) 

where x > 1. Define Z = lnx, multiple lnx to both sides, one has 

E E ^ - ' m = 1 - - E ^ ±Un)l (M) 

p r n 

which is valid for Z > 0. Note the analog with the classical trace formula [38T |. 

p r P 

with Z = t\[E. Thus one has formally 

oo oo 

tr£* = 1 + ^g n e-i» l = 1 -J2^ ±Un)l - (37) 

n n=l n=l 

This is very similar to the case for the motion on the surface of constant negative curvature, 
where the quantum eigen-momentum spectrum is {p n }- There tr£* = 1 + X^nLi e~^^ iPn ' 1 . Fur- 
thermore Biswas and Sinha have shown that the classical resonances are given by 7„ = | ± ip n [3f| 
by demonstrating that the peaks of the power spectrum of classical correlations indeed coincide 
with {p n }. 

For the RZ, the result of performing the auto-correlation of the smoothed density of states then 
leads to a self-dual "classical" spectrum {r;±it n }. So one might say that the classical resonances are 
7 n = \ ± it n . This interpretation of RZ as RP resonances was suggested recently 2) 3|- However 



here g n = —1 for the RZ while g n = 1 for the Riemann zeta function pole and RP resonances. 
This is related to the well-known sign problem for which a possible explanation has been offered 
by Connes 0. . 

We have clearly established a self-duality which then strongly suggests the interpretation of 
the spectrum {t n } as momentum or wave vector eigenvalues of a Helmholtz operator on a suitable 
domain. A similar description has been made in which the Riemann zeta function appears in the 
scattering matrix of the quantum scattering of a particle on a 2D surface of constant negative 



curvature 



2(], 2^|, so that {^tn} are the real parts of the poles of the scattering matrix. 



VI. CONCLUSION 



In this paper, we report the observation of new correlations among the RZ that presented as a 
conjecture supported by analytical arguments and numerical simulations. The spectral correlations 
possess three important properties: invariance, resurgence and prophecy. The prophecy can be used 
to progressively obtain new RZ from known ones. The last observation suggests an interesting new 
strategy for proving the Riemann hypothesis via a bootstrap approach starting with known zeros. 

We thank F.Y. Wu and J.V. Jose for discussions. This work is supported partially by NSF- 
PHY-0098801. 
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